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Abstract 



Following the work by Houghton, Reeve and Wallace about an alternative 
formulation of the n — > limit of the (n + 1) state Potts model in field theory 
for the large order behaviour of the perturbative expansion, we generalise their 
technique to all n by establishing an equivalence in perturbation theory order 
by order with another bosonic field theory. Restricting ourselves to a cubic 
interaction, we obtain an explicit expression (in terms of n) for the large order 
behaviour of the partition function. 



Resume 



Suivant le travail de Houghton, Reeve et Wallace concernant une formulation al- 
ternative de la limite n — > du modele de Potts en theorie des champs pour le 
comportement aux grands ordres, nous generalisons leur technique a n quelconque 
en etablissant une equivalence ordre par ordre en theorie de perturbation avec une 
autre thorie de champs de bosons. Nous restregnant a une interaction cubique, 
nous obtenons une expression explicite (en fonction de n) pour le comportement 
aux grands ordres de la fonction de partition. 



Introduction 



The Potts model[l[] is one of the most fascinating topics in statistical mechanics. It 
can be regarded as the generalisation of the Ising model to a generic number q of 
components. Originally this problem was proposed to Potts by his PhD Professor 
Dombf||], who suggested him to investigate the critical properties of a system of 
interacting spins confined in a plane, pointing in q equally spaced directions 

6, = — z = 0,l,...,g-l (1) 

with an hamiltonian depending only on the relative orientation of two neighbouring 
spins of the form 



w = - £ JCei.e.OoccosCei-ej) (2) 

{spins} 

Actually, the previous model bears the name of planar Potts model. Nevertheless, 
the most famous model is the one where two spins interact only when they are present 
in the same state, that is 

j(e 4 ,e,) oc s^, (3) 

The interaction can be re-written in such a way as to reflect a full-symmetry 
in a (q — l)-dimensional space by making use of q vectors (q — l)-dimensional e a 
pointing in the q symmetric directions of a hypertetraedron in (q — 1) dimensions. 
Indeed, the interaction (|3|) can be re-written as 

5(Q t , 0,) = -[l + (g-l) 6 e *.6 e ']. (4) 

The increasing attention to the Potts model was due to the fact that it has proven 
to be related to a large class of outstanding problems, especially when particular 
values of q, even non integer, are taken. For example, the percolation problem can 
be formulated in terms of q — 1 Potts model j|, [5|, |6|. Also Fortuin and Kasteleyn|7| 
showed that the Kirchhoff 's solution|| for an ensemble of resistors is strictly related 
to q = limit of the Potts model. Moreover, the critical behaviour of the dilute spin 
glasses is well understood in terms of the q — | limit [|, [10|] . Finally, Lubensky and 



Isaacon ]TTJ showed that interesting processes in gelation and vulcanization fall in 
the same class of universality of < q < 1 Potts model (For a complete review of 
the Potts model see F.Y.Wu|lj§). 
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A convenient formulation of Potts model as a field theory was first performed 



by Zia and Wallace ||13||, who investigated the critical properties in the framework 
of the renormalization group. Actually in the field theory literature it is customary 
to refer to (n + 1) Potts model (n + 1 = q in the previous notation). The sym- 
metry of the model allows a trilinear interaction which is the leading term in the 
renormalization group framework. Because of a pure cubic interaction, the high or- 
der behaviour of the perturbative expansion is clearly non-oscillating and therefore 
non Borel-sommable. Nevertheless, when particular limits are taken (e.g. n — > 
percolation problem) the large order behaviour is found to be oscillating and Borel- 
sommable |14]]) as obtained by Houghton, Reeve and Wallace (from now on referred 



as HRW). In Sec. I, we recall the hamiltonian of the (n + 1) Potts model and the 
alternative approach to n = limit used by Houghton, Reeve and Wallace. In Sec. 
II its generalisation to each value of n is showed. Sec. Ill will be devoted to the 
evaluation of the large order behaviour for such an alternative field theory and the 
results will be shown to agree with the previous ones. 

1 The fn+ 1)— state Potts model 



(5) 



We consider the [n + 1) state Potts model in d dimension, described by an 
Hamiltonian |13| of the form 

H[4>] = J d d x -(V0i)(V&) + + ^Pijk<f>i<t>j<f>k 

where sum over repeated indices is implied and runs from 1 to n. The system 
possesses a hypertetrahedric symmetry, provided that is of the form 

n+l 

p»k = £ (6) 

a=l 

where the e's satisfy 

n+l n+l n+l 

E e ? = £e«ej = (n+l)<% £ efef = (n + 1)^ - 1 (7) 

a=l a=l i=l 

Physical quantities can be extracted from the partition function (the discussion 
about the right contour in the path integral for a pure cubic theory can be found in 
[ jnp expanded in powers of g: 

„ oo 

Z n (g) = / V<f>e- H = £ Z K (n)g 2K (8) 

J K=0 
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We are interested in the asymptotic behaviour of Zx(n) as K goes to infinity. We 
can extract Zxin) from @ by the usual dispersion relation [[16], |i~7j 



Z K (n) 



-H[<f>i 



(9) 



where the closed integral is performed in the cut complex plane g 2 . In the large K 
limit, the previous integral can be evaluated by the steepest descent method in the 
space ((f>i,g) and the saddle point equations are 



1 

2 



- £ijk<Pj<Pk 



2K 



, , d d xe 
3! J ' g 



Equations ([U]) and (|TT1) can be decoupled easily by setting 



(10) 
(11) 



with the Mj's satisfying 



Equation flT3|) has n solutions u[\ with r = 1, 2, n, 



(12) 



(13) 



1=1 



where 



(14) 



,('') 



1 



(n + l)(n + 1 - 2r) 
Working in d = 6 dimension, the saddle point equations ([10|JTT| ) then read 



(15) 



V 2 </> c (x) = <j>l{x) 



(16) 



Y^UiUi J d & x(f>l{x) = -Kg 2 



(17) 



whose general solution is 



4> c (x) 



24A 2 



[\ 2 (x-x ) 2 + l] 2 

where A and x are constant parameters. From equation (O, 113) , we obtain 
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/ = 37T 3 ^rmnfe4 r M r) ) (19) 
and therefore ^ 

Z K(n) ~ K\ [ _ o n \_ (r) w | (20) 



3 • 2 8 7r 3 min{ r } ( u\ 'u\ 

For the (n+ 1) state Potts model (n > 1 integer) the minimum in ( p0|) is for r = 1, 
giving, as expected for a pure cubic theory, a malign behaviour for the coefficients 
Zk which do not oscillate [TBI, 15 1 



\2f m 1\2\K 



5n + l>-l , . 
Z «W~ Kl [ 3 .2Vn J (21) 

For n = (corresponding to percolation problem) the saddle point cannot be 
obtained by a naive analytic continuation of (|20| ) because evidently r = 1 is not the 
dominating saddle point. 

HRWfllifl found an alternative formulation for the percolation problem and they 



pointed out that in the n — > limit the Feynman rules are the same as those of the 
hamiltonian with two scalar fields 

i(V0) 2 -i(W) 2 + i(0 2 -^ 2 ) + 



H{M,g) = / d b x 



+ y9(<f> + ^) 3 



(22) 

provided that one keeps only (^—connected diagrams as those shown in figure 1. 
By means of path integrals techniques, they found the coefficients of the expansion 
of the two-point correlation function 



G(x, y)=Y,G K g 2K G K ~ X K\(-—) (23) 



which fortunately show an oscillating behaviour. 



2 HRW method for arbitrary n 

As mentioned in the introduction, it would be useful to generalise the previous treat- 
ment to each value of n. Our first task is to write down the equivalent hamiltonian 
in analogy with (p2|). In order to gain some insight, let us compute the two-point 
correlation function up to fourth order 
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(0e x tPijfc0j0i0fcPimn0«0m0nPafe C 0a060cPde/0d^e0/0cxt) (24) 

by performing all the possible Wick contractions. Because of the properties (|7|), we 
can get terms as 

n+1 

Y, e? [(n + l)S a ^ - l)[(n + 1)5 Q/3 - l][(n + 1)8^ - l]x 

a, /3, 7, 5=1 

x [(n + l)^" 5 - 1] [(n + l)cF 5 - 1] e 5 f (25) 

Inspired by the previous discussion, we try to interpret (n+l)5 af3 and (—1) as propa- 
gators of two different (possibly multi-component) scalar fields, and ifi respectively. 
The several terms we obtain by expanding the products in (EBp correspond in the 
second theory to all the different diagrams with the same topology as the original 
ones. As in the n — > case only diagrams with external legs connected by a 0's 
path do contribute. 

** Here Figure 1 ** 

In figure (1) we have drawn three of the possible graph built by and -0 propa- 
gators. By requiring that the diagrams in the two theories have the same numerical 
coefficients, we obtain the Feynman rules. 

Specifically, from the term involving only fields (figure (la)) we get 

n+1 

Y e? (n + if 8°^ 8^8^ 8^ P 6 e] = 

a, f), 7, 5=1 



n+1 

in+lfY^j (26) 

a=l 



In these case, the presence of 5's eliminate the sum over the four indices. If we 
associate a factor (n + 1) to each propagator, we obtain the same results without 
introducing any weight for the vertices. 

Let us now consider the second graph (see figure (lb)) 

n+1 

Y + lf8 ap 8 ps 8^8 aP {-l) 4 = 

«,/3,7,5 

n+1 

= (-l)(n + lfY (27) 

a=l 
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as we said, we associated a factor (—1) for the ip propagator and we do not need 
any factor for the vertex (f) 2 ip. 

The situation is quite different when we deal with graphs where we have a ifj 3 
vertex as in figure (lc) which derives from 

n+1 

(-1)(-1)(-1) E e?(n + l) 2 3*ty*ej = 

a,/3,7,5=l 

n+1 /n+1 \ 

= (-1)(^+1) 2 E E^/ • (28) 

7 =1 \a=l / 

We note that we are left with a sum over 7 which gives a factor (n+1). So a 
ip 3 vertex requires a factor of (n + 1). 

At higher order we also have non-zero diagrams with also <j) connected sub-graphs 
as shown in figure (2) 

** Here Figure 2 ** 

It can be rather simply proven that each 0— connected subgraph gives a factor 
of (n + 1). Summarising the Feynman rules for (n + 1) component Potts model are 

1) Only 0— connected graphs give a non zero contribution. 

2) Each <fi propagator gives a factor of (n + 1). 

3) Each ip propagator gives a factor of (—1). 

4) Each -?/> 3 vertex gives a factor of (n + 1). 

5) Each 0— connected subgraph gives a factor of (n + 1). 
These Feynman rules correspond to the following hamiltonian 



+ 



By following exactly the same procedure as HRW's in order to extract the 
)— connected diagrams, we get for the two-point correlation function 

(2) _ J V<P<P(x)<P(y) ex V j-H*(<P)]e nW * 



- / V^V(t) {l) V(p {2) (p {1 \x)(t) {2 \y) x 



cxp 



X 



-W{<j)^)-H*{4>W) 



e (l-n)W [J 2>0 e^+^w*] 
where TC*(<fi) is the reduced hamiltonian 



(30) 



= ((V0) 



1 1 /,_,, 2 , ,x I 1 



2n 



T ((V^) 2 + ^ 2 ) + 



and W* (■?/>) is defined by 



+ + # 



(31) 



exp[W* (?/>)] = / V<j>exp-H*((f>,ip) 



(32) 



3 High order behaviour 



As in the n — > case, it can be shown that only the second term on the left hand 
side of ( j3~0|) contributes to the high order behaviour of the perturbative expansion 
because of the particular form of the cubic interaction and we have to evaluate 



V^V^V^ (1) (x)0 (2) (y) — idg 2 -^- 

2m J q 2 l 



exp 



x 



g*K 

- W{(j) {2) ) + 2K\ogg 



(33) 



e (l-»)W [J x>ip e ( n+1 ) w *] 
for large Our goal can be achieved as usual by the steepest descent method. 
The equations for the saddle points are 



1 



71 + 1 

2 



V 2 ^ 



1 



9 (V> + & 



1,2 



1 1 <9W* 

— tVV + -^(^ + 0i) 2 +-g(^ + (j)2) 2 + (l-n)— — 
n+1 2 2 



1 

3! 



rf 6 x {(V + 0i) 3 + + 2 ) 3 } + (1 - ra)^ + — = 
1 J og g 



(34) 
(35) 
(36) 



It is simple to evaluate the partial derivative of W* with respect to g and ip by 
steepest descent method too and we obtain 



aw* 

dtjj 



— ^- V> - \g (0o + VO 2 ) l + 0(^ 2 ) 

71+1 2 
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and 



dg 



-~ (/^ + 0o) 3 ) l + 0(g 2 ) 



provided that <f) obeys to the following equation 



— ^rV 2 o = \g{<h + i>) 2 
n + 1 2 



After the re-scaling 



% = 0,1,2 



if; — ► — -ip 
9 



the saddle point equation are finally 



n + 1 



+ tf) i = 0,1,2 



V 2 ^ + -($! + + -($ 2 + 



^(*o + *) 2 



2if/ = J { j^ x {( $ i + + ( $ 2 + - (l-n)($o + ^) 3 } 



In order to solve the equations (0,0, |39|) , let us make the ansatz 



(37) 
(38) 
(39) 



= UiSJx) \1> = ud)Jx) 



(40) 



where <p c {x) is the solution of (|T6|). Therefore, we have reduced our differential 
equations to algebraic ones. 



———Mi = - [y + Mi 
(n + 1) 2 



■Ml 



1 (1-n) 

u 2 + 7 T-7T u 



(n + 1) (n + 1) (n + 1) 



3 • 2V 



6^3 



(z/ + Wl ) 3 + (i/ + n 2 ) 3 + (n - + «o) ; 



Solving the first equation gives two solution 



u 



(±) = - 



1 \ 1 mI 

±- -[l-2z;(n + l)] 2 



(n+1)/ (n + 1) 



Only four choices are possible 



(41) 
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1) 


If! 


= u (+) 


»2 


= «0 


= uH 


2) 


If! 


= uH 


«2 


= «0 


= u (+) 


3) 


Uq 


= ri (+) 


Wl 


= w 2 


= uH 


4) 


Uq 


= «H 


Ml 


= «2 


= «(+) 



From the 1) and 2) we obtain two equivalent solutions provided that 1 

2n 

z, = — n ^ 1 (42) 

(n + l)(n-l) 2 T y ! 

The choices 3) and 4) gives other two equivalent solution for v when n ^ 3 

4(1 -n) 

n^3 (43) 



(n + l)(n-3) 2 
The previous two value for v give two saddle points for g 2 



., 1 3-2 8 tt 3 n 



9 U = K^ (n + l) 2 (n-l) 2 ^ (44) 



and 



2 1 3-2 9 tt 3 (n-1) ,„ ,. cN 

^-5- (n + l)(n-3)^ ^ 3 (45) 
A this point, it is interesting to note that the saddle points g%\ and g%. are 
respectively the r = 1 and the r = 2 points of fll8l) . Anyway, Now, our final aim 
is to evaluate the large order behaviour of the perturbative expansion, that is the 
coefficients Zk in the expression (§). So a little care needs to understand which 
saddle point dominates at n fixed. We can obtain the two following behaviours 



7 (i) 

J K 



(n + l) 2 (n - l) 2 \ K 



3 • 2 8 tt 3 n 



n ^ 1 (46) 



and 



For n = (percolation) the second saddle point dominate over the first one and we 
recover the HRW results. Actually there exist other formulations for the percolation 



problem disagreeing with HRW ||18|| , but it can be shown Hl9| that the different 
result was due to an unappropriate choice of an integration contour. Nevertheless, 
the nature of the percolation behaviour still remains controversial because there is 
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another result [PD[] where the r = 1 saddle point seems to dominate. For n = 1 (Ising 
model) we cannot choose the coefficient Z$ and Z% is identically zero according 
to the fact that the trilinear interaction vanishes as n = 1. For n integer we have to 
take into account the ( pE6"| ) and this is in agreement with the calculations performed 
in section (|I|) 

Conclusions 

We performed the computation of the large order behaviour of the perturbative 
expansion of the (n + 1) state Potts model with a cubic interaction for a generic 
value of n, by mapping the original problem in a simpler one. We showed that 
our equivalent field theory has two saddle points depending on n and one has to 
choose the smallest one in order to extract the dominant term at large order. A 
natural question can arise looking at the different number of the saddle points in 
both the formulations. Indeed in Sec. I it has been shown that the Potts model, 
when studied by the usual techniques, has n saddle points characterised by an 
integer r, whereas our alternative formulation gives only two inequivalent saddle 
points. Indeed, the large order behaviour is related to the cut discontinuity of the 
correlation function for small value of the coupling constant, that is the imaginary 
part of the Green functions generated by the Taylor expansions around the saddle 
point. For n = 2, 3, ... the imaginary part of the correlation functions is dominated 
by the r = 1 saddle point. When an analytic continuation in n is performed and 
some particular limit are taken , the r = 2 saddle point becomes the leading term. 
All the other saddle points are sub-leading and they do not affect the large order 
behaviour. This generalisation can help to compute critical exponents for theory 
where n is not integer as those mentioned in the introduction eventually by adding 
higher order interactions. 

Acknowledgements 

It is a great pleasure to thank Giorgio Parisi for having suggested and supervised 
this work and G.F.Bonini for scientific discussions. 

References 

[1] R. B. Potts, Proc. Camb. Phil. Soc. 48, 106 (1952). 



12 



[2] E. Ising, Z. Phys. 21, 613 (1925). 
[3] C. Domb, J. Phys A 7, 1335 (1974). 

[4] P. W. Kasteleyn and C .M. Fortuin, J. Phys. Soc. Jap. (Suppl.) 26, 11 (1969). 

[5] M. R. Giri, M. J. Stephen and G. S. Grest, Phys. Rev. B 16, 4971 (1977). 

[6] H. Kunz and F. Y. Wu, J. Phys. C 11, LI (1978). 

[7] C. M. Fortuin and P. W. Kasteleyn, Physica 57, 536 (1972). 

[8] G. Kirchhoff, Ann. Phys. Chem 72, 497 (1847) 

[9] A. Aharony, J. Phys. C 11, L457 (1978). 
[10] A. Aharony and P. Pfeuty, J. Phys. C 12, L125 (1979). 
[11] T. C. Lubensky and J. Isaacon, Phys. Rev. Lett. 41, 829 (1978). 
[12] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 
[13] R. K. Zia and D. Wallace, J. Phys. A 8 1495 (1975). 

[14] A. Houghton, J. S. Reeve and D. J. Wallace, Phys.Rev. B 17, 2956 (1978). 
[15] A. J. McKane, Nucl. Phys. B 152 166 (1979). 

[16] H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and OPolymer 
Physics, (World Scientific 1990). 

[17] J. C. Le Guillou and J. Zinn Justin eds., Large Order Behaviour of Perturbation 
Theory, (North-Holland, 1990). 

[18] A.J.McKane, J. Phys. A 19 453 (1986) 

[19] S.Caracciolo and D.M.Carlucci, to be published 

[20] E.Brezin and A.J.McKane, J. Phys. A 20 2235 (1987) 



13 



List of captions 



Figure 1: Some diagrams built by means of ip (solid line) and ip (dotted line) propagators. 
Only the diagrams (a) and (6) are suitable because <p connected. 

Figure 2: High order diagrams with a <p connected subgraph. The filled region means 
any non-zero diagram. 
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